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ABSTRACT 

To examine the evidence for hierarchical evolution on mass scales of ~ 1O 13 -1O 14 M0, 
we apply a statistic that measures correlations between galaxy velocity and projected 
position (Dressier & Shectman 1988) to data for six poor groups of galaxies, HCG 42, 
HCG 62, NGC 533, NGC 2563, NGC 5129, and NGC 741. Each group has more than 
30 identified members (Zabludoff & Mulchaey 1998ab). The statistic is sensitive to 
clumps of galaxies on the sky whose mean velocity and velocity dispersion deviate from 
the kinematics of the group as a whole. The kinematics of galaxies within ~ O.l/i -1 
Mpc of the group center do not deviate from the global values, supporting our earlier 
claim that the group cores are close to virialization or virialized. We detect significant 
substructure (at > 99.9% confidence) in the two groups with the most confirmed 
members, HCG 62 and NGC 741, that is attributable mostly to a subgroup lying 
~ 0.3-0.4/i _1 Mpc outside of the core. We conclude that at least some poor groups, 
like rich clusters, are evolving via the accretion of smaller structures from the field. 
With larger poor group surveys, the incidence of such accretion and the distribution 
of subgroup masses are potential constraints of cosmological models on mass scales of 
< 1O 13 -1O 14 M and on physical scales of < 0.5h~ x Mpc. 
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1. Introduction 

The evolution of structure on different mass scales is one of the outstanding 
issues in cosmology. For example, although galaxy clusters of ~ 1O 15 M0 (including 
Virgo (Binggeli 1993; Bohringer et al. 1994), Coma (Mellier et al. 1988; Briel et 
al. 1992; White et al. 1993), and Abell 754 (Zabludoff k Zaritsky 1995)) are clearly 
evolving from the accretion of smaller groups, it is uncertain whether poor groups 
of ~ 1O 13 -1O 14 M also evolve hierarchically There is some indirect evidence that 
the evolution of poor groups is similar to rich clusters; the galaxies and hot gas in 
poor groups follow the same relationships found among the X-ray temperature, X-ray 
luminosity, and galaxy velocity dispersion for rich clusters (Mulchaey & Zabludoff 1998, 
hereafter MZ98). Historically, however, the number of known poor group members 
has been too small to examine individual groups for direct evidence of hierarchical 
evolution. 

Multi-object spectroscopy now makes it possible to obtain "cluster-size" samples of 
galaxies in poor groups and to identify substructure, if it exists, in the same manner 
as for rich clusters. Substructure in clusters was not detected until the number of 
spectroscopically-confirmed cluster members exceeded ~ 30-50 galaxies. Recent poor 
group surveys have reached this membership level (Zabludoff & Mulchaey 1998ab; 
hereafter ZM98a and b). The discovery of substructure in poor groups would provide 
new cosmological constraints by establishing that hierarchical evolution is occurring on 
mass scales of ~ 1O 13 -1O 14 M0 and on physical scales of ~ 0.5/i _1 Mpc. The detection 
of substructure would also support the picture in which at least some poor groups 
evolve as low-mass analogs to rich clusters. 

In this Letter, we describe the results from applying a substructure statistic 
(Dressier & Shectman 1988; hereafter DS88) to the six best sampled poor groups in 
ZM98ab, the most detailed spectroscopic survey of poor groups to date. 



2. The Group Sample 

In this analysis of substructure in poor groups, we consider six nearby systems 
(3800 < cz < 7000 km s" 1 ) drawn from our spectroscopic survey of 12 optically- 
identified groups (ZM98a). We choose the six groups, HCG 42, HCG 62, NGC 533, 
NGC 741, NGC 2563, and NGC 5129, because each has more than 30 spectroscopically- 
confirmed members (see ZM98a for the details of the membership algorithm). The 
group members range in absolute magnitude from the three or four brightest (M ^ M*) 
galaxies identified in past surveys to the large population of dwarfs (M & M* + 4) 
discovered in ZM98a. All six groups have X-ray halos extending out to radii of 
~ 200-300 kpc and X-ray temperatures of ~ 1 keV (MZ98). In this Letter, we use the 
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galaxy velocity data from ZM98a combined with new velocities for 51 additional group 
members (ZM98b). The number of known members within ~ 0.5/i _1 Mpc (~ the virial 
radius) of each group center ranges from 35 to 63 galaxies. 



3. Results and Discussion 

To search for substructure in the six poor groups, we use the method applied by 
DS88 to rich clusters. This test identifies a fixed number of nearest neighbors on the 
sky around each galaxy, calculates the local mean velocity and velocity dispersion 
of the subsample, and compares these values with the mean velocity and velocity 
dispersion of the entire group. The kinematic deviations of the subsamples from the 
global values are summed. This sum is larger for a group with a kinematically distinct 
subgroup than for a similar group without substructure. 

For each galaxy i, the deviation of its nearest projected neighbors from the 
kinematics of the group as a whole is defined as 5i = (n 1 / 2 / 'o r )[(vi oc — v) 2 + (cr r ; oc — 
^v) 2 ] 1 / 2 ) where v is the mean velocity for the group, ay is the group velocity dispersion, 
and n is the number of nearest neighbors (including the galaxy) used to determine the 
local mean velocity v\ oc and local velocity dispersion a r ^ oc . The total deviation for the 
group is defined as the sum of the local deviations, A = |<5j| for all i < N grp , the 
number of group members. As pointed out in DS88, the A statistic is similar to the 
X 2 statistic, except that the <Vs are not squared before summation in order to reduce 
the contributions of the largest, rarest deviations. If the galaxy velocity distribution of 
the group is close to Gaussian, and the local variations are only random fluctuations, 
A ~ ?V 

To calculate 5, we choose n = 11 (as in DS88). This choice allows robust 
determinations of v r> i oc and o r ,ioc- Silverman (1986) argues that using n ~ Ng^ 1 / 2 
nearest neighbors (= 6-8 for these groups) maximizes the sensitivity of such a test to 
small scale structures while reducing its sensitivity to fluctuations within the Poisson 
noise (also see Bird 1994b). To check the robustness of the n = 11 assumption, we 
compare the results below with those for n = 6. The conclusions drawn from the n = 6 
and n = 11 cases are the same. 

Calibration of the A statistic for each group is required because 1) the <5j's are 
not statstically independent and 2) the velocity distribution may not be intrisically 
Gaussian even if there are no subgroups (e.g., the group members may follow 
predominantly circular or radial orbits). We determine the significance of the observed 
A by comparing it with the results of 1000 Monte Carlo trials in which galaxy 
velocities are drawn randomly from the observed distribution and assigned to galaxy 
positions. This scrambling technique effectively destroys any substructure (DS88). If 
the probability is low that a group without substructure has a A value at least as large 
as that observed, then we consider the substructure detection significant. 
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In two groups, HCG 62 and NGC 741, the observed value of A is significant at 
the > 99.9% confidence level. Such high A values might arise from substructure, but 
also could result from smooth variations in the group's velocity field (e.g., rotation or 
velocity shear (Malumuth et al. 1992)) and/or from a dependence of ay on radius (Bird 
1994b). ZM98a show that a r is constant out to radii of ~ 0.5h~ l Mpc in the combined 
velocity dispersion profile for the sample groups. To determine whether substructure 
is in fact responsible for the high A values in HCG 62 and NGC 741, we examine the 
local deviation 5 for each group member. A concentration of large 8 values on the sky 
indicates a kinematically distinct subgroup. 

Figure 1 shows the projected spatial distribution of group members for each group 
(top panel). The second panel shows this distribution with the radii of the circles 
weighted by e 5 (as in DS88). Because each point is not statistically independent, a 
few very deviant galaxies can boost the S values for a large number of nearby points. 
Therefore, a visual comparison with the Monte Carlo simulations is required to assess 
the significance of any structures. The third panel shows the results of the Monte 
Carlo trial (out of 1000) with the largest A value, or greatest total deviation. The 
results of the trial with the median A value are in the bottom panel. 

The significance of the seven large clustered circles to the northeast of HCG 62 and 
five large clustered circles to the south of NGC 741 is high. In each case, the large S 
values show that a subgroup not in equilibrium with the global group potential is the 
principal source of the significant A value. Each of the two subgroups lies a projected 
distance of ~ 0.3-0.4/i _1 Mpc outside of the group core. 

On the other hand, the clustering of small circles within ~ O.l/i -1 Mpc of all of the 
group centers indicates that the core mean velocity and velocity dispersion are similar 
to the global values for the group. This result suggests that the group cores are close 
to virialization or virialized and is consistent with the conclusions from our earlier 
studies of group dynamics (ZM98a, MZ98). 



4. Conclusions 

Of the six best sampled poor groups of galaxies (each has more than 30 known 
members), two have significant substructure. In each case, the substructure detection 
is due mostly to a subgroup lying ~ 0.3-0.4/i _1 Mpc outside of the core. It is worth 
noting that we detect substructure in the two groups, HCG 62 and NGC 741, with the 
most confirmed members and also that our substructure statistic (from DS88) is not 
sensitive to subgroups that are well-superposed on the sky (Bird 1993, 1994a). These 
two points suggest that the fraction of poor groups with substructure is higher than 
the ~ 30% reported here. We conclude that poor groups, like rich clusters, evolve 
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hierarchically and that some poor groups are still accreting smaller structures today 0. 
In contrast, the cores of all of the poor groups (within ~ 0.1/i _1 Mpc) appear to be 
relaxed, an observation consistent with our earlier claims that the cores are close to 
virialization or virialized (ZM98a, MZ98). Deep spectroscopic surveys of groups in a 
large group sample would place limits on both the incidence of substructure and the 
distribution of subgroup masses, providing a new test of cosmological models on mass 
scales of ^ 10 13 -10 14 Mq and on physical scales of ^ 0.5/i _1 Mpc. 

We thank Dennis Zaritsky for his comments on the text. AIZ acknowledges support 
from the Carnegie and Dudley Observatories, the AAS, NSF grant AST-95-29259, and 
NASA grant HF-01087.01-96A. JM acknowledges support provided by NASA grant 
NAG 5-2831 and NAG 5-3529. 



Even the Local Group, which is less massive than these groups (~4x iO 12 M [Zaritsky et al. 1989] versus 
~ 10 13 -10 14 A/ o [ZM98a]), is forming from the collapse of the Milky Way and M 31 subgroups (Katm & Woltjer 
1959). There is also some evidence that AWM/MKW poor clusters, which are typically ~ 2-4 times more massive 
than these groups, have bound clumps at much larger radii \-2h~ 1 Mpc; Beers et al. 1995). 
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Figure Captions 

Figure 1: For each group, the top panel shows the distribution of group members on 
the sky. There are 35 galaxies in HCG 42, 63 in HCG 62, 36 in NGC 533, 41 in NGC 
741, 44 in NGC 2563, and 38 in NGC 5129. The dashed line in the NGC 741 and NGC 
5129 fields is the survey boundary (see ZM98a). The second panel is the same as the 
top panel, except that the radii of the circles are scaled by e 5 , where S is the deviation 
of the nearest neighbors from the kinematics of the group as a whole. The third panel 
shows the results of the Monte Carlo trial (out of 1000) with the largest value of the 
Dressler-Shectman (DS88) statistic A. The bottom panel is the trial with the median 
value of A. The scale bar below the bottom panel is 0.5/i -1 Mpc. 
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